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^ ! Abstract 

' In this work we study the particle production in time dependent 

(<~^ ■ periodic potential using the method of complex time WKB (CWKB) 

! approximation. In the inflationary cosmology at the end of inflation- 

^ I ary stage, the potential becomes time dependent as well as periodic. 

O . Reheating occurs due to particle production by the oscillating inflaton 

I , field. Using CWKB we obtain almost identical results on catastrophic 

^1 particle production as obtained by others. 

>< 1 INTRODUCTION 



The importance of particle production in expanding universe has been a mat- 
ter of interest and investigation since long times. Recently the catastrophic 
particle production has gained importance espcially in discussing the theory of 
reheating after inflation due to oscillating inflaton field. The literature on the 
subject, [|I], H, ^ I, 1^, H, H, |10|, |lT| though has contradictory statements, has 
been extensively reported in ref.[6,7]. The mechanism behind the catastrophic 
particle production is now called parametric resonance and stochastic reso- 
nance. A beautiful discussion of classical parametric resonance will be found 
in Landau and Lifshitz |T^. The recent attempts try to cast the equation of 
motion of the particle (that are produced) in the form of the Mathieu equation 
[p!3| starting with a basic model describing the inflaton field with a potential 
V{(j)) interacting with a scalar field x that are produced in the description. 
Earlier attempts to treat parametric resonance to explain reheating after in- 



flation were due to Dolgov and Kirilova 0, but in [15] the treatment was 
not rightly placed. 



From a few years back we are advocating |16, 17, 18, 19, EG, 21] a method to 



treat particle production in curved spacetime. The method is called complex 
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trajectory WKB (CWKB) approximation. The method is also apphed to the 



construction of wavefunction of the universe in quantum cosmology |2y, |2T 
with remarkable success and is very transparent from the standpoint of physical 
arguments. As mentioned already, due to differences in the growth rate factor 
even in the recent works [9,10,11], we take up the present work to study the 
reheating mechanism in inflationary cosmology through the method of CWKB 
and compare our results with the others. 

The basic principles and the mechanism inherent in CWKB are discussed 
in [13,14,15]. We discuss this in section 2 as a preparatory to the following 
sections. In section 3 we use the CWKB method to study particle production 
in the large amplitude region. This section will help us understand the ana- 
lytic theory of parametric resonance which is discussed in the section 4. We 
synopsize our findings in the concluding section. 

The basic idea that reheating occurs due to particle production by the os- 
cillating inflaton field was proposed by Linde |p. The particles so produced 
through interaction among themselves attain a state of thermal equilibrium 
with some temperature T. This process continues so long the scalar field 
transfers its all energy to the already produced elementary particles. The 
temperature at this stage is termed as reheating temperature, T^. A detailed 
discussion in this respect will be found in the beautiful papers H and ref- 
erences of previous works will also be found in that reference. 

Consider the decay of a scalar field (supposed to be inflaton) (p ^ XX- Iii 
the system if already there are many x particles with > 1, then the prob- 
ability is greatly enhanced due to Bose statistics. For fermion decay we will 
find that such enhancement would occur due to multiple refictions and will be 
explained in a separate publication. Because of periodicity emerging through 
the inflaton field (p the system evolves with explosive particle production. The 
temporal equation of x in a fiat Friedman background with scale factor a{t) is 
given by 

a f k'^ \ 
Xk + + {^ + ml{0)~^R+ = 0. (1) 

Here the time dependence of the background field and the scale factor is 
obtained from the evolution equations 

dV 

<P + ^HcP + — = 0, (3) 

where H = ^, V{(f)) is the effective potential of the scalar field. In a particular 
model, we can parametrize x equation in Minkowski space {a{t) = 1) as 

Xk + (k^ + g^cT^ + 2g^a^sinmt)xk = 0, (4) 
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where the periodicity = ^sinmt comes from the inflaton field. This 
equation is now of the form 

Xk + u;l{t)xk^O, (5) 

where the time dependent frequency is calit) — k"^ + g^o^ + Ig^asinmt. This 
periodicity is in the root for vigorous particle production and show parametric 
resonance for modes with certain values of k. The equation. (4) can be cast 
in Mathieu equation form 

Xfe + K - 1qcos2z)xk = 0, (6) 

with mt = 2z-n/2, = 4^^±^, q = • Here prime denotes difi^erenti- 

ation with respect to z. The properties of Mathieu equation shows that within 
the set of resonance bands of frequencies the modes grow as Xk ^xp{iJ,^^^z) so 
that it corresponds to exponential growth of occupation numbers of quantum 
fluctuations as nk{t) oc exp{2^^^\) and is interpreted as catastrophic particle 
production. Wc will now study the equation (6) to develop an analytic theory 
of parametric resonance. Before that we review the basic principles of CWKB. 



2 BASIC PRINCIPLES OF CWKB 

Consider the eqn.(5) and assume that u}{t) is nowhere zero for real t but when 
considered as a function complex t it has complex turning points given by 
1^(^1,2) = 0. According to the CWKB, the classical paths contributing in a 
complex semiclassical approximation joining two prescribed real points t' and 
t" (see Fig. 1) are composed of two parts coming from the contribution of the 
direct trajectory and the reflected trajectory. With 

Sitf,ti) = l\{t)dt, (7) 

the direct trajectory contribution is written as (we assume that t' » t") 

D.T=^=e'^^*"^*'\ (8) 

The reflected trajectory starts from t' and moving backward in time gets re- 
flected from the turning point ti and moving forward arrives at t" . This contri- 
bution is then multiplied by 0, 1, 2, 3, ... reflections between the the two complex 
turning points ti and ^2- This contribution is written as 
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Direct trajectory ,no reflection 
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K/ Reflected trajectory, one reflection 
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Figure 1: Trajectories with no reflection, one reflection and multiple reflections. 



Using 



iS{ti,t2)' 



H=0 



2n 



l + exp{2iS{ti,t2)y 



and taking t" — t we get with the replcement — > oo 

Xk{t, oo) ^ e'^(*'~) + iRe-'^^*'°°\ 



(10) 



(11) 



where we have neglected the WKB pre-exponential factor for convenience. The 
reflection amplitude is given by 



I g2iS(ti,t2)' 



(12) 



Using Feynman-Stuckleberg prescription and the boundary condition that 
there is no particle at i — > — oo i.e., 



Xk 



(13) 
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we identify 

2i5(ti,oo) I 

as the pair production amplitude with |i?cP + l^cP = 1- In (12) and (14) 
R (= full S-matrx element) is related to disconnected propagator and Rc (= 
connected S-matrix clement) is related to connected propagator [19]. The 
boundary conditions (11) and (13) arc known as scattering boundary condition 
and is applied in time dependent gauge. 

Hence for certain complex values of multiple reflection terms, (14) may 
exhibit poles for some parameter values if 

S{t,M)-{N+]^)T: (15) 

with N an integer. When the conditions (15) and uj{ti^2) = are not satisfied 
or the reality condition on z are violated, the physical-region pole becomes 
a resonance. This resonance occurs in parameter space and we will call it 
parametric resonance. The poles add up to give nonperturbative contributions 
enhancing Rc- Thus in CWKB we have a transparent idea of resonance particle 
production. Depending upon the parameter values in (15) (i.e., the values in 
Afe and q in eqn.(6)) we have resonance particle production in CWKB. In the 
next section we will discuss the particle production in large amplitude region. 



3 PARTICLE PRODUCTION IN THE LARGE 
AMPLITUDE REGION 

In the large amplitude region the parameter q is large (see eqn.(6)) and hence 
the resonance is broad. In that case we replace the periodic barrier by piece- 
wise inverted harmonic oscillators and within the interval |z| << 7r/4 we write 
cos{2z) ~ 1 — 2z'^. Let us try to understand this replacement. Considering (6) 
as a Schrodinger-like equation not in space but in time, in the potential problem 
qcos2z behaves as potential energy and Ak/2 as the total energy. For «\ 
q I, there is an infinite succession of large potential barriers separated by 
allowed region given by A}. — 2qcos 2z > 0. Instead of considering the barriers 
at once, we consider a single barrier and calculate the particle production 
which occurs at the turning points given by 

u;{z) = Ak- 2qcos2z = (16) 

These turning points occur at 

^ ^ ( (2n + l)7r ^ A ^^j^ ^ ^ ^2, ±3, ... (17) 
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for large q. Thereafter we repeat the process n times with n 'r^ zjix, n being 
a large integer. Now to calculate the particle production we are to determine 
particle and antiparticle states i.e., the left and right moving mode solutions 
with respect to which the in vacuum is defined. 

In CWKB these should be the adiabatic modes at large z. In other words 
the mode solution should satisfy the adiabaticity condition. As we have re- 
peated barriers we are to find out the adiabatic modes at finite z = z^ such 
that the modes effectively describe the vacuum at finite for large q. Let us 
consider the central barrier with — |(f + ^) < z < |(f — For large q 

the approximation gets better and the adiabatic condition is satisfied at finite 
z = zq ior \ z tt/A but far away from the turning point z = |(7r — ^) 
(considering a single barrier). In such a case the above replacement is a good 
approximation. As we are considering adiabatic modes far away from the 
turning points (basically these are the WKB modes) , the reflection and trans- 
mission coefficients or the Bogolubov coefficients will be constant for a barrier. 
The maximum of potential barriers occurs at vr, 2tt, Stt, ... and hence the reso- 
nance occurs for a broad range of values of k with k"^ /m? = — 2q. Now the 
resonance width varies as q^ with 1—1,2, 3, ... . So for narrow resonance, the 
first band with I — 1 will be important with Ak — 1 ± q. Let us consider the 
broad resonance. In this section we consider only the production of particles 
at each barrier as if there is no particles in the previous barrier and also ne- 
glect the enhancement at the zeroes of the infiaton field where non adiabatic 
transitions occur . Introducing cos2z ~ 1 — 2z^ and with 

z' = (4g)VS, A = V^e, e = ^ - 1, (18) 

the Mathieu equation (6) reduces to the form 

^ + {X + z")^^0. (19) 

The left and right moving modes having the asymptotic behaviour of unit flux 
[10] are given as, 

XL,R = (2^;2)(-i^^^^/2)exp(Ti^^2) 

which can be obtained from the solutions of (19) written in terms of parabolic 
cylinder functions as 

XLiz) (X L>_(,^,+i)/2(e-/42gV4^) 

with 

Xr(z) ^Xl{z)- 
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Despite the finite range zq << it/ 4, the large q in the argument of D-function 
makes it possible to use the asymptotic form of these D-functions to use as 
mode solution as if we are considering the WKB adiabatic modes. In CWKB 
we do not require the exact mode solutions to evaluate the particle production 
amplitude; however we are to guarantee that the in vacuum is rightly specified. 
To evaluate R as in (12) we are to evaluate the complex turning points. The 
complex turning points of (19) are at z' = Using CWKB expression of 

reflection amplitude, we find after the evaluation of the integrals 



2 e"^^ 

\Rr.\ 



1 + e-^^ ' 



(20) 



When we consider full Mathieu equation, the turning points will be repeated 
and the solution will be of the form 

V'(2) ~ iR^Piz) = e^' P{z) = e"'"l^lp(z), (21) 

where n ~ 2;/7r is the number of barriers crossed during the time z. Hence 

1 



TT 



ln\Rc 



1 exp{-n^e/2) ^^2) 
TT 1 + exp{—Ti y/qe) 



If q is large, we get 



A - -v/^e/2 (23) 

This result exactly coincides with the result obtained by Fujisaki et. al.||ll||. 
This result amply reflects the usefulness of the complex turning points in the 
CWKB method. The construction of CWKB wavefunction as ~ |-R|" is how- 
ever a crude estimate. In obtaining (23), the multiple reflections between the 
complex turning points have been ignored and the periodic aspect of the po- 
tential has been introduced in a qualitative way through (21). However, a 
general treatment, taking the effect of repeated barriers, will now be taken 
within the framework of CWKB. 

It should be pointed out that while calculating the reflection potential as 
I R I" we basically assume that (i) there is no particles produced in the previous 
barrier , (ii) and neglect the non adiabatic changes at the zeroes of the infiaton 
field and the phases that the field acquire before entering a potential barrier 
due to scattering. 

It should be pointed out that in the above discussion we have obtained the 
reflection and transmission coefficients as constants when we consider adiabatic 
evolution far away from the turning points that lie around the vicinity of the 
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zeroes of the inflaton field. To take into account the non adiabatic changes 
that occur around the zeroes of the inflaton field, we will expand uJ^{t) around 
the zeroes of the inflaton field. For the purpose we would bring u;'^{t) term 
in the form A + Bsin?mt and then expand around the zeroes t = tj of the 
inflaton field. In this case also we can safely assume the adiabatic evolution 
in the region away from the zeroes t = tj of the inflaton field. We can then 
take the coefficients of the CWKB adiabatic modes exp{-^iuj{i)) as constants 
until we face the next barrier where the coefficients will change. We take up 
all these aspects in the next section. 

4 ANALYTIC THEOTY OF PARAMETRIC 
RESONANCE 




Figure 2: The CWKB trajectories in periodic potential 

Our objective in this section is to develop an analytic theory of parametric 
resonance particle production using the method of CWKB. We consider a 
Schrodinger equation with time dependent potential 

^ + [A;^ - Vm = 0, (24) 

where the periodicity in the potential is given by V{t + T) = V{t), where T 
is the period. In the previous section we have not considered the effect of 
periodicity explicitly. When V{t) is periodic the situation is quite interesting. 
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We now have two types of reflection points, one is complex and other is real in 
time. For the latter we will use tunneling boundary condition akin to particle 
production in space dependent gauge |2^. Suppose V{t) = Vocos2t with 
period T = tt. The solution of (24) will then behave as exp [±(A;^ — Vq)^''^] as 
if at these points the solutions undergo non-adiabatic changes. This implies 
that between two such reflection points, as clarified in the previous section, 
the solution undergoes adiabatic evolution where the WKB approximation 
can safely be used. We now call the points tj where V(tj) = Vq as 'reflection 
points' and the complex points where k"^ — V{t) = as the 'turning points'. 
The inflaton which generates this V{t) is zero at these points, (see (4)). 
There is another aspect of periodicity. In a particular barrier, between tj and 
tj+i, there will be some particles already produced with respect to the previous 
barrier and hence the boundary condition of 'no particle state' at t — * — oo 
will not be satisfied. This fact is to be taken into account while considering 
particle production in a given barrier. 
Let us define 

u;{t) = [e - V{t)Y/^ (25) 

and ^ 

S{tf,ti) = j ^ uj{t)dt (26) 

as the action for the path that goes from to to t. Suppose we have two com- 
plex turning points ti and t2 given by cij(ti 2) = 0. We adopt the boundary 
conditions such that at i — *• —00 we have the transmitted wave 

^/^(t^-oo) =Tee*^(*'*o), (27) 

and at t — s> +00 

ilj{t +00) = e^^(*'*o) + i?ee-*'^(*'*o^ (28) 

Here Rc and are respectively the reflection and the transmission amplitude. 
In CWKB the expression of the reflection amplitude is given in (14) and Tc 
is determined from the relation |-RcP + l^cP = 1- Equation. (27) and (28) 
have the interpretation that at t ^ —00 we have no particle but at t ^ +00 
we have pair production in the out vacuum with respect to in vacuum. The 
same problem can be evaluated in terms of Bogolubov mode decomposition 
technique with 

ij{t ^ +00) = «^e'^(*'*°) + /3^e-^^(*'*o), (29) 

in which the the Bogolubov coefficients are given in terms of transmission and 
reflection coefficients as [20] 

M = J^2^ 1/5^1' = W (30) 

I c| I c| 

In the above treatment the particle and antiparticle states are defined at t ^ 
±00 with the corresponding vacua respectively as out and in vacuum. We now 
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proceed toward the construction of CWKB eigenfunction. Let tj and tj^i be 
the two points where V{tj) = and t is a point such that tj<t< tj+i where 
we want to calculate the eigenfunction. We choose the rightmoving and left 
moving waves as follows (see fig. 2). 

Rightmoving wave from to to t — exp [—iS{t, to)] , 
Leftmoving wave from to' to t — exp [iS{t, to')] , 

where S{ti,tf) is given in eqn.(26). For the rightmoving wave the transmission 

and the reflection coefficients arc (Tk,Rk) whereas for the left moving wave 
we take the coefficients as {T^,Rl). We now avoid the subscript c on i? and 
T. Let in the region I, tj_i < t < tj, the CWKB solution before tunneling is 
(represented by (1) and (!') in Fig. 2) 



The Bogolubov coefficients in a given barrier are now supposed to be constants. 
Between tj < t < tj+i, in the region II we take similarly 

^l+\t) = ^e-^^(*'*o) + i^fc e+^^(*'*o), (32) 

where the coefficients and f^l'^^ are constants for tj < t < tj+i. Hence- 
forth we will not write the WKB pre-exponential factor for convenience. The 
constuction of ^;^^^(t), according to CWKB is now shown in the fig. 2. In the 
region tj < t < tj^i the rightmoving part consists of two parts. The part (1) 
represented by (xlexp[—iS{t,to)] after transmission gets multiplied by l/T^ 
because of transmission at tj and using (30) we find 

Qi-j 

ai exp [-iS{t.to)] exp [-iS{t, to)] . (33) 

Tk 

There is also a contribution to the rightmoving wave coming from the region 
III where t > tj, as shown in the fig. 2, giving the part (3lexp[+iS{t,to)] in the 
region t < tj represented by (!') in the figure 2. The amplitude of the left- 
moving part at tj is Pi exp [+iS{tj,to)]. This amplitude part when continued 
in the region t > tj becomes Pi/T^ which again being refiected at tj becomes 
(5{RI/T^ so that in the region t > tj 

Pi exp [iS{t, to)] ^ 4? exp [+iS{tj, to)] exp {-iS{t, tj)). (34) 
-'fc 

Writing S{t,tj) = Sit,t,) + S(t.iJ,o) - S{tj,to) = S{t,to) - S{tj,to), we get 
from (34) for the rightmoving part as 



i'k,RMP 



I Pi 2iS{tj,to) 

Tk n 



^-iS{t,to)_ (35) 
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To get (35) we can also start from the region III with (3l'^'^exp{iS{t,t'Q)) that 
now gets multiphed by Rl/T^ because of transmission at tj+i and reflection 
at tj and then use the continuity condition 

We will get the same result as (35). Comparing (35) with the first term in (32) 
we get 



— + 



(36) 



Now the left moving part has two contributions given by the trajectory (3) 
and (4) of fig. 2. The trajectory (3) coming from t'^ gives the left moving part 
jSj^'^ exp {iSit^t'o)) in the region t < tj+i. Now we use the relation above to 
convert j3l^'^ in terms of This part then in t > tj becomes 



1^1 iS{t,to) 



(37) 



Another contribution comes from exp {—iS{t, to)) which on transmission at 
tj and reflection at tj^i gets modified to 



—R f.-iS{tj+i,to)+is(t,tj+i) 



(3^ 



Now, iS{t,tj+i) = iS{t,tj+i)+iS{tj+i,to)-iS{tj+i,to) = iS{t,to)-iS{tj+i,to) 
so that the left moving part becomes 



■ip' 



i+i 

k,LMP 



Pi 



+ 



Q^fc Rk 2iS(t,+i,to) 



M5(t,to) 



(39) 



Hence comparing (39) with the second term in (32) we get 



Pi. 



n 



k ^-2iS{t,+i,to) 



(40) 



This result exactly coincides with Kofman, Linde and Starobinsky [0]. We 
have in the expression of jSl'^^ the phase term as S{tj^i, to) instead of S'(tj, to)- 
It should be pointed out that we have not taken repeated reflections between 
the turning points tj and t^+i which will automatically be introduced when we 
calculate Rk and Tk through the technique of CWKB. 

Now we are to determine the transmission and reflection coefficients using 
CWKB in the region tj < t < tj+i. For the purpose we need to specify 
V(t). We take V(t) = —g'^(^'^siin?{mt). Actually this type of terms arise 
in inflationary scenario from an effective potential ^(0) = ^0^, and the 
interaction term — |(7^0^x^- Instead of (4) we now have 



Xk + (fc' + g^^^ sin^{mt))xk = 



(41) 
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In the vicinity of tj equation (41) is transformed to 

^ + + g'<^>'m\t - t,r) X, = 0. (42) 
With T ^ K{t - tj), \ = , kl ^ g^m, (38) reduces to 

^ + (A^ + r^)- (43) 
Using the results in (18) and (28) we find 

R, = J^-^ , (44) 
Vl + e^ 

In CWKB the phase (pk is unknown. However for the particular problem it 
can be calculated knowing the solutions of (39) in terms of parabolic cylinder 
functions and we get 

0fc = ar^rf^— j. (46) 

We now calculate the number density of outgoing particles as n;^"*"^ = l/^fc^^P 
and find after simple algebra 

n^^+i = e-'^^' + (l + 2e-'^^')ni 

-2e-5^Vl + e-'^^V"'fc(l + ^i) ^^^^iou (47) 

where the phase Q\oi — 2^^"*"^ — 0fc + CL'f'gf^i — argal- In obtaining (47) we have 

used \ai\ — = 1 andjo;^! = y^I + nj,). Here the arbitrary point to is now 
taken as zero and 

ei+^ = S{tj+u 0) = r^\{k'' + g''^^sin\mt))]^/^dt. (48) 

t/ 

For narrow resonance i.e.,g = << 1 and with Ak — ^ + 2q we will find 

S(tj+i, tj) ~ ^/A'k7^ = nn, (49) 

taking into consideration the effect of periodicity i.e., at the points tj the phase 
in (51) must be an integral multiple of tt. We have also neglected the terms hke 
0{q/ Ak) while evaluating the right hand side of (48) through eUiptic integrals. 
Thus 

e^^^r = (-i)Vs (50) 
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and our result then exactly coincides with the results given in [6,7,8]. Apart 
from this minor difference all the results arrived in [6,7,8] will also be the same 
in our case. Our CWKB treatment in this section and in the previous section 
(eqns. (22) and (23)) simply reflects the region where both the results are 
valid. 

Let us briefly mention the important results obtained from analytical anal- 
ysis. 

(i) In CWKB the resonance particle production occurs due to rotation of cur- 
rents at the 'turning points' and 'reflection points'. The former is independent 
of time and we call it spontaneous particle creation and the latter is due to 
multiple reflections at the points tj and we call it induced particle creation. 

(ii) Due to spontaneous particle creation the number of produced particles 

always increases whereas the induced creation may result in a destructive in- 
terference between tj and tj^i when sinOj remains positive and varies. This 
solely occurs due to the time dependence contained in 9j and multiple reflec- 
tions between tj and tj+i. 



(iii) Actually the spontaneous particle creation determines the resonance 
structure. A large ttA^ will suppress the the effect of particle creation. Hence 
we must have 

where Ak = ^ + 2q , q = This gives an estimate of resonance width as 

<kl^ gm^/n. (52) 

(iv) Now exp{—TiX'^) oc exp{—l/g) and is non-analytic at g = 0. Thus (45) 
manifests the non-perturbative nature of CWKB resonance effects. 

(v) In the broad resonance regime q » 1. We take ^{t) — constant; it 

makes 9k and (pk independent. It can be shown that in that case the analytical 
solution nicely matches with the numerical solution of Mathieu equation. The 
details will be found in [6,7,8]. 

Let us consider an interesting interpretation of CWKB reflection coefficient. 
From the expression (12), has the following heuristic interpretation. Since 

2 ^ ^ (53) 
' ' n+l(n + l)' ^ ^ 

and l/(n -|- 1) is the probability that in vacuum evolves into the out vacuum, 
we say that out of n + 1 particles in the \out > vacuum, n particles get reflected 
so that gives the absolute probability of one pair production from the in 
vacuum as if n paricles remain in the \in > the vacuum. Hence (1 + n)^|i?p 
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gives the number of produced pairs in the \out > vacuum. Now in CWKB 
exp{—7rX^) is the number of particles produced in any barrier. The number nl. 
particles already are in the in vacuum, so we do not require the first factor in 
(53). Thus 

\R\' = -T^ni ^ nl (54) 
n;^ + 1 

Hence 

(l + e-^^'yni (55) 

gives as the number of particles produced in the out vacuum while nl. remaining 
in the in vacuum. We recall that p3i 



|^|2 _ g-/n(l+JV(fc))|^^|2 _ g-2/m£.„(fcm|^j2_ ^gg) 

When no particles are produced N{k) = 0, we get |-R| = \Rc\ and ImCeff = 
leading to \in >= \out > . In that case we have two possibilities: (i) the ini- 
tial vacuum contains no particles and hence \R\ = \Rc\ = 0, (ii) the initial 
vacuum contains some particles (squeezed \in > state) then from (12), (1 + 
exp{2iS{h,t2))y\R\'^ = {l + expi2iS{ti,t2))Y\Rc\'^ = (l + exp(2i^(ti, ta))) VP 
gives the number of particles in the out vacuum due to presence of n particles 
in the squeezed \in > vacuum. Hence the number of particles ni^^ in the 
barrier between tj < t < tj+i due to nl particles in the |m > vacuum is 



nl = particles produced in out vacuum for zero particle in state 
+ {l + expi2^S{t^,t2))Y\nl\'). 

Hence, 

4+1 = e-^^' + (1 + e'^^fnl (57) 

This result coincides with the first two terms of (47) when exj9(— vrA^) << 
1. U nl = i.e., zero particle in the in vacuum we get the standard result 
nk = exp(— 7r)A^ for any barrier. The skeptical readers may avoid this heuristic 
arguments. 

Next, consider the large occupation number situation in which n^ » 1. 
In that case we write (45) as 

nl=nlexp{2nfii), (58) 

where 

1^1 = —In (1 + 2e'"^' - 2smCe-"^'/Vl + e"-^'). (59) 
2% 

After a number of infiaton oscillations we write the number of x particles 
occupation number n^ as 

nfc(t) = ie'"^^ ^'^ ~ e^"^/* (60) 
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so that 



^ 47r2 



/ dkk^e^"""-'. (61) 



Wc now evaluate the integral by steepest descent method with /jLmax = at 
k = km and using n'l ~ 2/1/ Sk'^, we finally get 



Stt^ ^/tt/i mt 



~ Z' (62) 



The above result comes out the same as that obtained in [7] . 



5 CONCLUSION 

In this work we find that the complex trajectory WKB approximation nicely 

reproduces all the essential aspects of particle production due to inflaton os- 
cillation. The reason of difference with other works [10, 11] is due to not 
taking the multiple reflections (a non-perturbative quantum corrections) in 
their approach. Except for a phase factor the coincidence of our results with 
that of [6,7,8] leads credence to the CWKB approach. In this work we have 
concentrated the treatment to the Minskowski spacetime only to understand 
the effect of periodic potential on particle production. With a{t) ^ 1, the 
qualitative conclusion will not change remarkably. In curved spacetime uJ^{t) 
will be modified as 

u}'^{t) = + g^^'^sin^mt + Aa; 

where 

Ao; = - |(a/a)2 - |(a/a) - 

is usually very small in the region of inflaton oscillation where « rn? 
and can be neglected. In this case we are to replace x by Xk{t) = a^^^Xfe(t). 
In the broad resonance regime the effect of a{t) might result in the number 
of produced particles at some time depending upon the parameter values and 
the FRW scale factor. The details in this respect will be reported elsewhere. 
Actually in our approach we do not require much to depend on Mathieu equa- 
tion, everything can be settled out from the CWKB framework. The stochastic 
resonance (in which the number of particles may decrease at some moments), 
backscattering and fermion production can also be taken in our frame work. 
In our previous works [15, 16] we have discussed fermion production without 
repeated reflections between the 'reflection points'. We like to take all these 
aspects in a future publication. 
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